Topological Magnetic Insulators with Corundum Structure 
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Topological insulators are new states of quantum matter in which surface states residing in the bulk 
insulating gap are protected by time-reversal symmetry. When a proper kind of antiferromagnetic 
long range order is established in a topological insulator, the system supports axionic excitations. 
In this paper, we study theoretically the electronic states in a transition metal oxide of corundum 
structure, in which both spin-orbit interaction and electron-electron interaction play crucial roles. 
A tight-binding model analysis predicts that materials with this structure can be strong topological 
insulator. Because of the electron correlation, an antiferromagnetic order may develop, giving rise 
to a topological magnetic insulator phase with axionic excitations. 

PACS numbers: 71.70.Ej, 75.30.Kz, 75.80.+q, 73.20.-r 
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The discovery of time reversal invariant topological in- 
sulator has attracted great attention in condensed mat- 
ter physics [3-0|. With time-reversal symmetry broken 
on the surface, the electromagnetic response of three di- 
mensional (3D) insulators are described by the topologi- 
cal 9 term of the form So — ^jrln I d 3 xdtE, ■ B together 
with the ordinary Maxwell terms, where E and B are 
the conventional electromagnetic field insides the insu- 
lator, a = e 2 /he is the fine structure constant, and 
is the dimensionless pseudoscalar parameter describing 
the insulator, which refers to "axion" field in axion elec- 
trodynamics For a system without boundary, all 
the physical quantities are invariant if 9 is shifted by 
integer multiple of 2tt. Therefor all time reversal invari- 
ant insulator fall into two distinct classes described by 
either 9 = (trivial insulator) or 9 = it (topological 
insulator) (l~0| . Such a universal value of 9 = tt in topo- 
logical insulators leads to magneto-electric effect with an 
universal coefficient, which has several unique experimen- 
tal consequences such as a topological contribution to the 
Faraday rotation or Kerr rotatio n II lHl3j . and the image 
monopole induced by an electron [14(. 9 has an explicitly 
microscopic expression of the momentum space Chern- 
Simons form which depends on the band structure of the 
insulator [111 
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where Af v {k.) = — i^^d / dki\u u ) is the momentum space 
non-abelian gauge field, with \uy), \u v ) referring to the 
Bloch wavefunction of occupied bands. 

If strong electron correlation exists in a topological in- 
sulator, a long-range antiferromagnetic(AFM) order can 
be established under low enough temperature. Since 
the AFM order breaks time-reversal symmetry sponta- 
neously, 9 can deviate from 7r, and also becomes a dy- 
namical field which has fluctuations associated with some 
spin collective modes. The spin collective mode induc- 



ing fluctuations of are thus coupled to the photons by 
#E ■ B term, which means they become "axions" in the 
term used in high energy physics. Such a nonconventional 
antiferromagnetic insulator supporting axion excitations 
is proposed as topological magnetic insulator (TMI) [l5| . 
Due to its coupling to photons, the axion field hybridizes 
with photons, leading to axion polariton, with a polariton 
gap tunable by an external magnetic field. Thus such a 
material can be used as a novel type of optical modulator 
to control the transmission of light through the material. 

To realize the TMI phase, we need both the nontriv- 
ial topology of the electron bands and strong electron 
correlation. The materials with electrons in 4d or 5d- 
orbital can have both strong spin-orbit coupling (SOC) 
and strong interaction, which is ideal for this purpose. 
Recently, models for topological insulators with strong 
electron correlation have been proposed [16|-|l9j] , also first 
principle calculations show topological phases exist in 
thallium-based III-V-VI2 ternary chalcogenides [20I 21 j as 



well as ternary heusler compounds which contain the rare 
earth element Ln, where additional properties ranging 
from superconductivity to magnetism and heavy-fermion 
behavior can be realized [22L [23j. In this Letter, we study 
theoretically the transition metal oxide ABO3 of corun- 
dum structure with B and A standing for some transition 
metals such as Fe, Ti, Ru, Rh, Ir, Os, etc [24J. A tight 
binding model is obtained by using point group symme- 
try of this structure, from which we find a topological 
magnetic insulator phase with certain SOC strength and 
electron-electron interaction. 

The corundum structure is shown in Fig. [lja). Each 
transition metal atom is surrounded by oxygen octahe- 
dron, and the d orbitals are split by the octahedral crys- 
talline field into doublet e g (x 2 — y 2 , 3z 2 — r 2 ) and triplet 
t 2g (xy,yz, zx) orbitals (See Fig. QJc)). We will neglect 
small distortion of the oxygen octahedra which may lead 
to minor corrections to electronic structure [25j. The 
energy of ti g stays lower with respect to e g , because 
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FIG. 1: (color online) (a) Corundum crystal structure with 
three primitive lattice vector denoted as ti,2,3. Each tran- 
sition metal ion M (M=Ir, Os, etc.) (green large circles) is 
surrounded by oxygen octahedron (red small circles). Each 
unit cell has four M atoms, denoted as 1,2,3,4. The space 
group is R3c. (b) Brillouin zone for corundum structure 
with space group R3c. The four inequivalent time-reversal 
invariant points are r(000), F(iv00), L(tvtttv) and Z(OOvr). (c) 
Schematic crystal field splitting of 5d level in corundum struc- 
ture. We are interested in the half filling case with effective 
angular momentum j e ff = 1/2. 



the latter point towards the negatively charged oxygens. 
The SOC is effective in t 2g orbitals and negligible in e g 
orbitals. Including the SOC, t 2g splits into total an- 
gular momentum j e g = 3/2 and j e ff = 1/2. We fo- 
cus on those materials where the Fermi level lies com- 
pletely in the j ff = 1/2 sub-bands. For example, the 
ions Ir 4+ , Os 3+ , Ru 3+ etc with five <i-electrons satisfy 
this requirement^. 

To obtain the electron dynamics in this system, we 
start by a symmetry analysis to the corundum structure. 
The space group of this structure is D 3d (R3c) with four 
atoms in each unit cell. It has a trigonal axis (three- 
fold rotation symmetry C3) defined by z axis, a binary 
axis (two- fold rotation symmetry C 2 ), denned by y axis, 
and inversion symmetry with the inversion center at the 
middle of the two neighbor transition metal atoms. The 
primitive lattice vectors £1,2,3 and primitive unit cells are 
shown in Fig. QJa) , where each unit cell consists of four 
transition metal atoms denoted as 1,2,3,4. Since the O 
p-level e p are far away from the fermi level, we can con- 
sider a model describing only d-electrons, with the hop- 
ping mediated by the oxygen p-orbitals. The model is 
generally written as 



Hq — 



[4 



tij dj 



E 



[4*ij d J 



h.c 



(2) 



where {ij) and ((ij)) denote the nearest-neighbor (NN) 
and next-nearest-neighbor (NNN) sites, respectively, and 
the hopping terms tij and t\j are in general 2x2 ma- 
trices. The form of the parameters Uj,^ can be sim- 



plified by symmetry considerations. Due to space lim- 
itation, we will only present the result of the symme- 
try analysis. The NN transfer integral tij are real and 
spin independent, with two independent parameters, the 
intra-plane hopping t and the inter-plane hopping t±. 
t = (pdn) 2 [(ppa) + 3(ppir)]/3(e d - e p ) 2 [13], where (pdn), 
(ppcr), and (ppir) are Slater-Koster parameters between 
pd and pp, respectively [27j |. The contribution of the 
order of (pdn) 2 /{td — e p ) cancel out in the honeycomb 
lattice, in sharp contrast to S^IrCU with the perovskite 
lattice [2^, (2t| . The NNN transfer integrals are spin de- 
pendent, and it is essential for the realization of the topo- 
logical insulator phase. For intra-plane in A plane, 1 — > 1 
hopping can be written as 

i' xl = it'^a ■ f u + ^11, (3) 

fn is a unit vector fu cx txx + l/s/22, txx is the hopping 
link. £22 = *ii d ue ^0 inversion symmetry. While in B 
plane, i' 33 = e^'lije^ 2 "' , i' u = 
due to C2 symmetry. For intra-plane (A — > B) } 
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fx3 is a unit vector r*i3 cx tx3 — az, tx3 is the hopping 
link, a is some parameter which depend on materials 
and cannot be determined purely by symmetry, below we 
choose a — l/y/2 which has almost the same amplitude 



as intra-plane. i' 24 = n3i t'x 
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Explicitly, for the intra-plane 1 — > 1, 2 — > 2 hopping 
are Xx,yx> z x, and 3 — > 3, 4 — > 4 are x%, y 2 , z 2 denoting in 
Fig. Ha). 

In summary, the transfer integrals are real and spin 
independent for NN links, while complex and spin de- 
pendent for NNN links. The accurate hopping parame- 
ters varies in different materials. As an example, in the 
following we will use the transfer integrals of Ir oxide 
introduced in Ref. 17]. One can always define all the 
parameters in the unit of in-plane nearest neighbor hop- 
ping t, which leads to t = 1, t^ = 0.33, tx\\ = —0.1, 
t± = y, t 2 ± = 0.5?/, t' 2X = Xt 2 ± = OAy. Here A is 
the SOC strength which determines the ratio of spin- 
dependent hopping and spin-independent hopping. For 
Ir oxide we have A = 0.8. All the inter-plane hopping 
matrix elements are rescaled by a factor y which incorpo- 
rates the anisotropy between intra-plane and inter-plane 
directions. The energy dispersion for y = 0.3 (dashed 
line) and y = 0.55 (solid line) are shown in Fig. [2a), 
which shows that the system at half filling is an insulator 
in both case. Due to inversion symmetry, all the energy 
bands are doubly degenerate. 

In three-dimensional topological band insulators, four 
independent Z 2 topological invariants can be defined [30I — 
32]. For inversion symmetric systems, all the topolog- 
ical invariants can be simply determined by the par- 
ity of the wave- functions at the 8 time-reversal invari- 
ant momenta (TRIM) in the Brillouin zone Q. De- 
note Gi, G2, G3 are the three basis vectors of the 
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FIG. 2: (color online) (a) 3D energy band dispersion of tight- 
binding model for corundum structure with A = 0.8, and 
y = 0.3(dashed line), y — 0.55(solid line), (b) The change 
of energy levels at F point (7r00) versus y. A band crossing 
occurs at y = 0.4245. The system changes from trivial insula- 
tor to topological insulator, (c) & (d) 2D band structure for 
a slab with 001 surface for the parameters A = 0.8, y = 0.3 
in (c) and y = 0.55 in (d). The red curves in (d) stands for 
surface states. The inset shows the surface Brillouin zone. 



reciprocal lattice, then the 8 TRIM's are denned by 
ki = (fciGi + k 2 G 2 + fc3G3)/27r with ^1,^2,^3 = or 
7T. For each TRIM k^, one can define a Z 2 quantity Si 
as the multiplication of the parity of all occupied bands 

k = rLeocc £ s ' wn " n £ s tne P ari ty °f s_ th band. It should 
be noticed that a Kramers pair of bands are only counted 
once, otherwise Si would always be even. The four Z 2 in- 
variants (uq] U1V2V3) can be determined by Si in which 
i>q = Y\. i= i Si is the strong topological invariant which is 
stable upon disorder, and responsible for the topologi- 
cal magneto-electric effect 11| . More discussions on the 
other three "weak topological invariants" can be found 
in Ref.[3]. For y = 0.3, we find 8 = +1 at all TRIM 
at half filling. On the contrary, for y = 0.55 we find 
5 = -1 at the three F points (see Fig. Q](b)) and 8 = +1 
at all other TRIM's. Consequently, y = 0.55 phase is 
a strong topological insulator(STI) with the topological 
character (1; 000), and y = 0.3 is a trivial insulator with 
character (0;000). From this result we see that a band 
inversion]!] occurs at F points upon the change of y. In 
Fig. [2] (b) we show the energy at F point versus y, from 
which one can see clearly a level crossing at y ~ 0.42. 
The topological invariants can be calculated for all val- 
ues of anisotropy parameter y and spin-orbit coupling 
parameter A, which leads to the phase diagram shown in 
Fig. [3] One can see that the topological nontrivial band 
structure can be realized at large y (i.e., small anisotropy) 
even for infinitesimal spin-orbital coupling. However, one 
should notice that for some parameters the band struc- 
ture is actually a semi-metal (similar to Sb), which has 
a direct gap but does not have in-direct gap. We also 
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FIG. 3: (color online) The phase diagram of the system with 
two variables: the out of plane hopping parameter y and SOC 
parameter A. The green and gray regions stand for topolog- 
ical nontrivial and trivial phases, respectively. Point A and 
B correspond to the parameters used in Fig. HJc) and (d), 
respectively. 



solve the Hamiltonian ^ in a slab geometry with two 
001 surfaces to study explicitly the topological surface 
states. Fig. (2fc)&(d) shows the 2D energy dispersion of 
the two systems shown in Fig. [2fa). In addition to the 
bulk states, for y — 0.55 there are surface states with 
three Dirac cones at M points of the surface BZ, while 
no surface state is found for y — 0.3, in consistency with 
the bulk topological invariants. 

To get better understanding of the physical proper- 
ties of this system, a low energy effective model can be 
obtained by expanding the Hamiltonian around the F 
points. Around each F point, the effective model is 4 x 4 
which describes two Kramers pairs of low lying bands and 
has Dirac-like form. In the following, we will denote the 
momentum by its coordinate in the basis of reciprocal lat- 
tice, i.e., k = (fcid + k 2 G 2 + k 3 G 3 ) /2tt. The F points 
are given by (it, 0,0), (0, tt, 0) and (tt,tt, 0). Around the 
point (tt, 0, 0) the Hamiltonian has the following form: 

5 

H cS (tt00) = e (q)I 4x4 + d a (q)T a , (5) 

a=l 

Here the Dirac T-matrices are defined as T a = (t x ® 
&x , t x ® Cy, T y ® 1 , i~z ® 1 , t x ® o z ) where r, and <7j 
(i = x, y, z) denote the Pauli matrices in the space 
of orbital and spin, respectively. q = k (tt00), 

d»(q) - Ei=i A aA°?< for a = M. 3 . 5 * d M) = m + 

T,i=i,2,s B i9i, and e (q) = C + Ei=i,2,3 D ^f- For 
A = 0.8, around the topological phase transition point 
we have 

/ 0.14 -0.12 0.37 -0.34' 
A? = -0.47 0.06 -0.13 0.09 
\ 0.014 0.038 0.015 0.055 

Bi = (0.625,0.32,0.24), 

A = (0.375,0.04,0.04) 



4 



C = 0.064. The mass parameter M depends on y as 
M ss — 23y + 9.76 which changes sign at y ~ 0.42 and 
leads to the topological phase transition. The effective 
Hamiltonian around the other two F points at (07r0) and 
(7T7r0) can be obtained by C3 rotation. 

Now we study the effect of electron correlation. The 
leading term in the interaction Hamiltonian is the onsite 
Hubbard repulsion for the j c g = 1/2 orbitals 
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Magnetic ordering in this system can be studied in mean- 
field approximation. For simplicity, we will only consider 
the order parameters that do not break translation sym- 
metry. The mean-field calculation predicts a SDW phase 
above a critical U, as shown in the phase diagram in Fig. 
21 The spin moments of this SDW phase lie in the hon- 
eycomb plane, which are ordered antiferromagnetically 
within each honeycomb layer and non-collinear between 
the two neighboring honeycomb layers, as shown in Fig. 

El 

Such a SDW order breaks both time-reversal sym- 
metry T and inversion symmetry V spontaneously, but 
preserves the combination VT ■ Thus the magneto- 
electric coefficient 9 in the 6*E • B term may deviate 
from the time-reversal invariant values and tt in the 
SDW phase, and certain spin wave fluctuations are cou- 
pled to photon as axions[l5j. To estimate the value of 
9 in the SDW phase, we can study the effect of the 
SDW order in the effective model (0. To the lead- 
ing order, the T, V breaking but TV preserving per- 
turbation to the effective model must have the form 
i5ff e ff(q) = J2 a =i,2.3,5 m a^a, where m a depends linearly 
on the spin moment (S s ) obtained from mean-field the- 
ory. Thus the perturbed Hamiltonian can still be writ- 
ten as H(q) = £)o=i d a (q)T a with d a (q)(a = 1, 2, 3, 5) = 
Y^i= x y z Afh + m a and di(q) unchanged. For this effec- 
tive model, 9 has an explicit expression (IBj 
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^did^dydAdl (7) 



where i,j,k,l = 1,2,3,5, and |d| = y^ =1 ^. Since 
the main contribution to 9 comes from the region close 
to Dirac points, 9 can be approximated by the sum of 
9s calculated separately for each Dirac point using the 
effective model. The numerical results of 9 is shown in 
Fig. 11(b). 
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FIG. 4: (Color Online) (a) The phase diagram with the out 
of plane hopping y and onsite repulsion U as parameters. 
The phase on the right is the topological magnetic insulator 
(TMI) which carries the dynamic axion. (b) The value of 9 
(see text) along the blue line in the phase diagram, (c) The 
SDW order pattern. The purple arrows represent the spin in 
the honeycomb layer A and the green arrows represent the 
spin in the adjacent honeycomb layer B. Other possible spin 
configurations can be obtained by six-fold rotations of this 
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